Y  -4D-A110  606  WASM1N6T0N  UN IV  SEATTLE  OCFT  Of  CHEMISTRY  r/9  7/H 

|  VIBRATIONAL  ENERSY  TRANSfER  ANO  NYROCYSIS  Of  NITRONE THANE  BY  TH~ETCtU> 

JAN  62  *  YUAN*  B  S  RA8IN0VITCH*  R  TOSA  N0001A-7B-C-0690 

'UNCLASSIFIED  NL 


1 


BJEILE  COPV 


& 


LEVEI/ 

Vibrational  Energy  Transfer  and  Pyrolysis  of  Nitromethane 
by  the  Variable  Encounter  Method 


by  W.  Juan,  B.  S.  Rabinovitch  and  R.  Tosa 


Department  of  Chemistry  B6-10,  University  of  Washington 
Seattle,  Washington  98195 


Technical  Report  No.  NR092-549-TR24 
Contract  N00014-75-C-0690,  NR-092-549 


/ 

January  15,  1982 


Prepared  for  Publication  in 
Journal  of  Physical  Chemistry 


OFFICE  OF  NAVAL  RESEARCH 
Department  of  the  Navy 
Code  473 
800  N.  Quincy 
Arlington,  VA  22217 


Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  '  ^  ^ 

the  United  States  Government.  This  document  has  been  approved  for  --  '  '  *  / 

public  release;  its  distribution  is  unlimited. 


f 


< 


Unclassified 

SECURITY  Cl  ASSlFICATiON  Jt  iHI'j  PAoL  />«(*  Fnlffwiij 


REPORT  DOCUMENTATION  PAGE 


1  REPORT NUMbfcR 

H 

NR092-549-TR24 


4  title  »u 

Vibrational  Energy  Transfer  and  Pyrolysis  of 
Nitromethane  by  the  Variable  Encounter  Method 


HKAIj  INS  I  KIK  TIONS 
HKFOKF.  COMF1  I-  TINT.  FORM 


RECIPIENT'S  CATALOG  NUMBER 


5  TYPE  OF  REPOR1  »  PERIOD  COVERED 


Technical 


(  PERFORMING  ORG  REPORT  NUMBER 


W.  Yuan,  B.  S.  Rabinovitch,  and  R.  Tosa 


9  PERFORMING  ORGANIZATION  N»M£  AND  ADDRESS 

Professor  B.  S.  Rabinovitch 
Department  of  Chemistry  BG-10 
University  of  Washington 

if Is .  WA _  98195 _ 

II  CONTROLLING  OFFICE  NAME  AnO  ADDRESS 

Office  of  Naval  Research,  Code  473 
Department  of  the  Navy 
800  N.  Quincy 

1  imrinn  \I  A  '7991 


iNu  AuEnCY  NAME  A  AODUE  S  Sf  I  f  dtllorapt  from  Controlling  Oltlco) 


6  CONTRACT  or  grant  NUMBERfs) 

N00014-7  5-C-0690 
NR092-549 


10  PROGRAM  ELEMENT  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 


12  REPORT  DAI  fc 

January  15,  1982 


1  3-  NUMBER  OP  PAGES 
20 


15  StCURITY  CLASS.  thlm  tmport) 

Unclassified 


15«.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


16  OliTRiauTlOH  STATtWtHT  thl  •>  K»po»t) 

This  document  has  been  approved  for  public  release; 
its  distribution  is  unlimited. 


17  DISTRIBUTION  statement  (ot  ihm  •barrtrr  «nl*r«J  in  Block  iO,  It  dUlmrant  from  /import) 

Prepared  for  publication  in  J.  Phys.  Chem. 


19  KEY  WORDS  (Coni  lima  «>»»  rmvmram  iM*  It  nacaaaary  and  Idantlly  by  block  numbar) 

Energy  Transfer  Surfaces 

Gases  Unimolecular  Reaction 

High  Temperature  Vibrational  Relaxation 

Nitromethane 


20  ABSTRACT  fConflnut  o»i  ravaraa  aldm  It  nacaaaary  and  Idantify  by  block  numbor) 

The  pyrolysis  of  nitromethane  has  been  studied  by  the  Variable  Encounter 
Method  (VEM)  at  temperatures  from  816  K  to  1092  K  with  two  reactors  of  differing 
geometry  having  fused  silica  surfaces.  The  probability  of  reaction  per  collisioi 
with  the  reactor  surface  was  measured.  The  down  energy  transition  jump  size, 
<AE'>,  was  determined.  It  decreased  with  increasing  wall  temperature.  A  com¬ 
parison  is  made  of  <AE’>  with  previous  results  reported  to  date  for  other  sub¬ 
strate  molecules.  Nitromethane  is  one  of  the  more  efficient  energy  transfer 
agents.  However,  in  addition  to  the  size  (vibrational  eigenstate  density)  and 
the  polarity  of  the  molecules,  the  nature  of  the  hot  surface  seems  also  to  play 


DD  ,  j, 


COITION  OF  I  NOV  t*  It  OBtOLCTE 
S/N  0102  LF  014  B601 


Unclassified 

iFCllBiTV  Cl  AtilFir  ATirtM  r\W  Tw»A  PA  r.W  rtotm  »*»•*•*#> 


Vibrational  Energy  Transfer  and  Pyrolysis  of  Nitromethane 
by  the  Variable  Encounter  Method 


t 

by  W.  Yuan,  B.  S.  Rabinovitch  and  R.  Tosa 

Department  of  Chemistry  BG-10,  University  of  Washington 
Seattle,  Washington  98195 


(  Abstract 

-I 

The  pyrolysis  of  nitromethane  has  been  studied  by  the  Variable 
Encounter  Method  (VEM)  at  temperatures  from  816  K  to  1092  K  with 
two  reactors  of  differing  geometry  having  fused  silica  surfaces. 

The  probability  of  reaction  per  collision  with  the  reactor  surface 
was  measured.  The  down  energy  transition  jump  size,  <AE’>,  was 
determined.  It  decreased  with  increasing  wall  temperature.  A  com¬ 
parison  is  made  of  <AE‘>  with  previous  results  reported  to  date  for 
other  substrate  molecules.  Nitromethane  is  one  of  the  more  efficient 
energy  transfer  agents.  However,  in  addition  to  the  size  (vibrational 
eigenstate  density)  and  the  polarity  of  the  molecules,  the  "at  e  of 
the  hot  surface  seems  also  to  play  a  role. 
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Introduction 

Nitromethane  is  a  widely  investigated  species  because  of  its  importance 
in  a  number  of  phenomena  including  the  mechanism  of  gas  phase  nitration  of 
hydrocarbons,  the  mechanism  of  formation  of  photochemical  smog,  and  the 
chemistry  of  propellants.  A  large  amount  of  work^’^  has  been  devoted  to  the 
study  of  the  kinetics  and  mechanism  of  the  thermal  decomposition  of  nitro¬ 
methane.  However,  virtually  no  attention  has  been  paid  to  the  study  of  vibra¬ 
tional  energy  transfer  involving  this  species.  Of  course,  coliisional  acti¬ 
vation  and  deactivation  is  a  fundamental  physical  process  of  great  ubiquity. 
Since  nitromethane  is  a  highly  polar  compound,  it  is  also  of  interest  to  study 

the  behavior  of  this  molecule  by  the  VEM  method  which  has  been  applied  recently 

7  8 

to  a  variety  of  molecules,  principally  non-polar  hydrocarbons.  ’  By  this 
technique,  vibrational  relaxation  in  the  transient  region  may  be  studied, 
in  principle,  on  a  collision-by-collision  basis. 

In  this  method,  a  molecule  equilibrated  at  some  low  temperature 
in  a  reservoir  flask  is  allowed  to  enter  a  reactor  of  variable  dimensions, 
say  cylindrical ,  heated  to  reaction  temperature  (called  an  encounter). 

Two  events  are  possible:  the  molecule  may  escape  from 
the  finger  after  a  sequence  of  collisions  with  the  hot  wall;  or  the  series  of 
vibrational  relaxation  collisions  with  the  wall  may  be  terminated  by  a  reaction 
event.  Only  subjects  for  which  the  decomposition  reaction  is  a  homogeneous 
process  have  been  chosen  for  study.  By  varying  the  geometry  of  the  finger, 
the  average  residence  time  of  the  molecule  in  the  hot  reaction  zone  can  be 
varied;  and  the  average  number  m  of  successive  collisions  that 
occur  before  the  molecule  exits  from  the  finger  (and  re-equilibrates  with  the 
cold  wall)  can  therefore  be  varied.  The  average  reaction  probability  per 
encounter  can  be  measured  experimentally  for  different  fingers,  and  a  suitable 


energy  transfer  model  may  be  estimated  by  fitting  the  data.  The  collision 
events  themselves  may  be  simulated  by  a  Monte  Carlo  calculation  using  the 
known  reactor  geometry; and  a  distribution  function  f(n,m)  may  be  found  for 
the  number  of  sequential  collisions  n  that  occur  per  encounter  with  a 
reactor  of  given  m  (m  =  h). 

In  this  paper  we  describe  the  study  of  nitromethane  in  two  silica 
reactors  at  temperatures  between  800  K  and  1100  K. 
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Experimental 

Nitromethane  (Aldrich  Chemical  Company,  99%,  spectrophotometric  grade) 
was  used  without  purification.  Gas  chromatographic  analysis  showed  that  it 
contained  less  than  0.5%  nitroethane. 

Q 

As  was  described  previously,  the  apparatus  consisted 
of  a  1-1  iter  fused  quartz  spherical  reservoir  flask  which 
had  two  quartz  cylindrical  finger  reactors  blown  on  to  the  surface.  These 
fingers  were  of  cylindrical  geometry  having  length/radius  ratios  of  'v  1.7  and 
'v  6.4,  corresponding  to  mean  numbers  of  sequential  collisions,  m  =  5.0  and 
14.4,  respectively.  Each  finger  could  be  separately  heated  in  a  stainless 
steel  furnace  block  which  was  maintained  by  ceramic  clamshell  Kanthal  electric 
heaters.  The  temperature  of  each  finger  was  measured  with  five  chromel-alumel 
thermocouples  cemented  separately  to  the  outside  of  the  reactor  finger.  Typical 
reactor  finger  temperatures  were  -v  800  K  to  ^  1100  K;  flask  reservoir  tempera¬ 
tures  were  400  K  to  480  K.  The  temperature  variation  during  the  run  time  was 
±  1 °C  and  along  the  finger  length  was  ±  10°C.  The  temperature  of  the  reactor 
was  estimated  by  simple  averaging.  The  effective  temperature  of  the  cold  flask 
was  suitably  weighted  by  surface  area,  but  calculation  showed  that  the  effect 

on  the  theoretical  calculation  due  to  a  change  of  the  flask  temperature  by  as 
was  negligible. 

much  as  50°Ka  Before  a  run,  the  reactor  and  flask  reservoir  were  evacuated  to 
approximately  3-5xl0~6  torr.  Reaction  run  pressures  were  in  the  range  7-8x10” 
torr.  Reaction  time  varied  from  several  minutes  to  ten  hours,  depending  on  the 
reaction  temperature.  The  reaction  percentage  varied  between  2%  to  40%.  At 
the  end  of  a  run,  the  reaction  mixture  was  expanded  from  the  reactor  into  a 
liquid  nitrogen  trap. 

Separation  of  the  reaction  products  was  made  at  100°C  on  a  1-m  column  of 
3-mm  copper  tubing  packed  with  Porapak  Q.  The  detector  was  FID. 
Acetonitrile,  which  is  inert  under  these  conditions. 
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was  used  as  an  internal  standard  against  which  the  disappearance  of  substrate 

could  be  accurately  measured.  No  attempt  was  made  to  monitor  all  of  the 

1-5 

product  species, which  for  this  system  are  very  numerous. 

Surface  Seasoning 

The  surface  condition  of  the  hot  finger  had  a  significant  influence  on 
the  observed  rate  constant.  If  the  surface  were  seasoned  with  reactant  itself 
as  has  been  customary  in  all  previous  work  -  a  technique  that  has  led  to  mini¬ 
mal  and  reproducible  rates  of  reaction  -  the  rate  constant  became  progressively 
higher  although  eventually  leveling  off.  This  suggested  the  possible  presence 
of  surface  catalysis.  It  was  found  that  oxidizing  the  surface  at  1100  K  with 
30-40  torr  pure  oxygen  for  several  hours  and  then  evacuating  the  reactor  to 
3  -  5xl0-6  torr  gave  reproducible  results  and  minimal  rates  of  reaction.  The 
oxidation  procedure  was  carried  out  at  each  temperature  studied. 


Results  and  Discussion 


The  thermal  decomposition  reaction  obeys  the  first-order  law.  This  is 
illustrated  in  Fig.  1.  The  apparent  rate  constants  are  sunmarized  in  Table  I. 
Each  tabulation  is  the  average  of  from  three  to  four  determinations  from 
separate  runs.  As  Fig.  2  shows,  the  Arrhenius  relation  (which  no  longer 
has  a  simple  meaning  here)  was  satisfied  for  this  temperature  range.  The 
apparent  activation  energies  calculated  from  Fig.  2  are  38.8  kcal  mole”1,  for 
the  m  =  5  finger,  and  42.4  kcal  mole'1,  for  the  m  =  14.4  finger.  The  latter 
value,  being  closer  to  the  low  pressure  steady  state  value,  should,  of  course, 
be  the  higher  one. 

There  exist  several  different  experimental  sets  of  high-pressure  Arrhenius 

1-6 

parameters  for  this  reaction.  We  chose  to  use  the  measurement  of  Glanzer  and 
5  -1 

Troe,  E,  =  58.5  kcal  mole  ,  which  is  quite  close  to  the  value  calculated  by 

a 

Benson  and  O'Neal6  from  thermodynamic  quantities  and  mechanistic  considerations. 
Since  this  activation  energy  is  the  highest  of  all  those  reported,  we  are 
inclined  to  believe  that  this  choice  is  optimal.  The  GT  and  BO  values  for 
the  frequency  factor  are  peculiarly  different,  but  these  quantities  have  no 
bearing  on  our  results;  only  Ea  enters,  i.e.,  the  critical  threshold  E  , 
since  our  data  were  taken  in  the  second-order  region. 

The  vibrational  spectrum  of  nitromethane  has  been  assigned  9  (see  App.  I). 
The  value  of  Eq  for  nitromethane  was  calculated  to  be  54.4  kcal  mole'1,  with 
use  of  the  average  temperature  1150  K  over  the  range  used  by  GT  and  the  activa¬ 
ted  complex  model  of  App  I. 

The  average  probability  of  decomposition  per  collision  with  the  hot  wall 
in  a  reactor  of  given  m-value,  Pc(m),  was  calculated  from  the  apparent  first- 
order  rate  constants,  the  known  reactor  dimensions  and  the  reservoir  flask 
temperature.  Since  our  rates  were  calculated  from  nitromethane  disappearance, 
it  is  necessary  to  consider  the  possible  extent  of  chain  or  sensitized  character 
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of  nitromethane  pyrolysis.  According  to  Perche's^  simplified  mechanism  used 
for  computer  simulation  kinetics,  the  important  secondary  radical  reactions 
which  cause  significant  decomposition  of  nitromethane  are 

ch3no2  +  ch3  — ►  ch4  +  ch2no2 

and 

ch3no2  +  ch3o  ch3oh  +  ch2no2 

From  the  detailed  product  construction  given  by  these  authors  at  700  K,  the 

additional  decomposition  fraction  of  the  parent  is  estimated  at  approximately 

3Q%.  Since  these  free  radical  reactions  have  low  critical  thresholds,  the 

te  lerature  increase  in  our  system  to  the  range  800  K-1000  K  favors  the  uni- 

molecular  decomposition  relative  to  these  abstraction  processes.  Moreover, 

in  our  system  the  total  pressure  is  very  much  lower  than  in  the  usual  pyrolytic 

conditions  so  that  it  is  expected  that  the  importance  of  second  order  reactions 

are  even  further  reduced.  In  fact,  our  rate  constants  were  unaltered  within 

-4 

experimental  error  when  the  reaction  pressure  was  increased  from  3x10  torr 
-3 

to  3x10  torr.  This  implies  that  in  our  system  the  observed  rate  constant  is 

virtually  independent  of  the  role  played  by  bimolecular  secondary  reactions 

and  we  have  adopted  the  observed  rate  values  as  the  correct  ones.  The  F>c(m) 

values  are  collected  in  Table  II  and  Fig.  3  for  m  =  5  and  14.4. 

The  computer  simulation  of  the  encounter  activation  process  has  been 
7  15 

described  elsewhere.  ’  Two  models  have  been  used  to  characterize  the  proba¬ 
bility,  p. .  of  a  down  energy  transition,  AE  ,  by  the  molecule  from  energy  E. 

■  J  J 

to  energy  E.. : 

Model  FG  (flat  gaussian,  i.e.,  AE  independent  of  Ej, 

Pij  =  A1  exp(-(AE-AEmp)2/2o2)  ; 

Model  BE  (exponentially  weighted  Boltzmann), 


Pi j  =  A2N(E.)  exp  (E./RT)  exp  (-AE/<AE>)  . 


The  correlation  table  used  in  the  computer  simulation  is  given  in  Table  III. 

The  correlation  coefficient  C.  is  defined  as  the  probability  that  molecules 

J 

which  have  already  experienced  j  collisions  with  the  hot  wall  will  have  at 

oo  oo 

least  the  (j+l)th  collision,  i.e.,  C.  =  l  n./  £  n-,  where  n.  is  the  sum  of 

J  i=j+l  1  i=j  1  1 

the  molecules  in  the  distribution  vector  left  after  i  collisions.  The 
correlation  coefficients  sumnarize  the  varying  probability  of  a  molecule 
remaining  in  the  reactor  after  each  wall  collision.  The  parameters  of  each 
particular  model  distribution  were  adjusted  to  provide  theoretical  agreement 
with  the  experimental  curve  for  the  m  =  5  finger;  the  parameters  so  chosen 
were  then  used  to  calculate  a  theoretical  curve  for  the  m  =  14.4  finger. 

Fit  is  made  to  the  m  =  5  curve  because  it  provides  greater  sensitivity  to 
model  detail  than  does  the  m  =  14.4  curve  which  lies  closer  to  the  steady 
state  value.  The  graphs  are  shown  in  Fig.  3.  It  can  be  seen  that  the  theo¬ 
retical  curves  are  close  to  the  experimental  value  for  m  *  14.4,  showing  good 
concordance  between  the  two  sets  of  experimental  data.  The  BE  model  fits  the 
data  somewhat  better  than  the  FG  model  but  there  is  little  to  choose  between 
them.  A  set  of  parameters  for  the  FG  model  that  exactly  fit  the  experimental 
m  =  14.4  data  may  be  calculated  and  give  somewhat  higher  values  of  <AE’>. 
Calculated  steady  state  curves  (m  =  °°)  as  well  as  the  strong  collider  curve  are 
also  given  in  Fig.  3,  and  it  appears  that  the  hot  wall  will  approach  the  strong 

collider  behavior  when  the  temperature  falls  substantial ly  below  800  K.  Once 

11-13 

again,  the  trend  observed  in  previous  VEM  studies  of  increasing  strength 
of  collisions  with  decreasing  temperature  is  borne  out.  (However,  the  lowest 
temperature  value  for  m  =  14.4  appears  somewhat  too  high.)  The  effective  aver¬ 
age  energy  of  down  transitions,  <AE’>  for  models  BE  and  FG  together  with  their 
parameters  are  summarized  in  Table  IV. 
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Comparison  with  the  previous  results  by  the  VEM  method  is  made  in 
Table  V.  A  rough  rule  that  the  down  jump  energy  transfer  coefficiency  per 
collision  of  the  substrate,  declines  with  increasing  size  (i.e.,  increasing 
vibrational  eigenstate  density)  of  the  molecule  seemed  evident  from  earlier 

Q 

studies  with  various  hydrocarbon  reactants.  Also,  it  was  shown  that  the  polar 
molecule,  iodopropane  (y  =  2.0  D) ,  seems  more  efficient  in  energy  transfer 
than  non-polar  ones,  due  to  the  longer  sticking  time  that  gives  rise  to  more 
complete  accommodation.  The  measured  efficiency  for  highly  polar  nitromethane 
(y  =  3.6  D),  while  one  of  the  more  efficient  in  Table  IV,  seems  a  little  lower 
than  expected  (which,  incidentally,  also  belies  any  significant  role  for  sur¬ 
face  catalysis).  However,  we  note  that  in  the  present  work  an  oxidized  surface 
was  used  instead  of  the  surface  seasoned  by  the  reactant  itself,  as  in  our 
previous  studies,  so  that  we  may  be  seeing  an  effect  of  different  energy  trans¬ 
fer  efficiency  of  the  surface.  Indeed,  as  mentioned  above  under  Experimental, 
the  rate  (not  reported)  on  a  seasoned  surface  was  larger;  calculated  values 
of  the  energy  transfer  on  such  a  surface  were  found  and  gave  magnitudes  that 
even  exceed  those  for  iodopropane. 

Clearly,  apart  from  the  polarity  or  size  of  molecules,  the  nature  of  the 
surface  must  play  a  significant  role  in  molecule-surface  relaxation  process. 

We  will  attempt  to  clarify  this  feature  in  future  studies. 
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Table  I.  Apparent  Rate  Constants 


m_  T  ( K ) 


10^  k  (sec~^ ) 


5.0  870  4.1 

5.1 

4.2 

987  59 

70 
64 
84 

1069  290 

290 
320 
320 

14.4  816  12.5 

10.1 

8.4 

12.2 

907  139 

143 
176 
133 

1013  2000 

1760 
1780 
1470 

1092  8200 

8000 
7800 


for  Decomposition 


Av.  4.5  ±  0. 6a 


Av.  69  ±  8 


Av.  305  ±  12 


Av.  10.8  ±  1.4 


Av.  148  ±  14 


Av.  1750  ±  100 


Av.  8000  ±  200 


a)  Standard  deviation  of  the  mean 
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Table  III.  Correlation  Coefficients  C 

n 


n/m 


5 


14.4 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


0.750 

0.759 

0.788 

0.833 

0.800 

0.889 

0.809 

0.907 

0.815 

0.921 

0.819 

0.930 

0.822 

0.937 

0.824 

0.942 

0.824 

| 

0.946 

i 

0.950 

i 

0.952 

0.954 

0.956 

0.957 

0.958 

0.958 

i 

i 

i 
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Table  IV.  Energy  Transfer  Parameters  for  Nitromethane 


T(K) 

816 

907 

1013 

1092 

FG 

<AE  > 
mp 

2440 

1960 

1760 

1700 

A 

> 

m 

V 

2700 

2170 

1950 

1880 

BE 

A 

> 

m 

V 

911 

943 

981 

1090 

U 

A 

L±J 

< 

V 

4970 

4100 

3360 

3355 

o 


a)  Energy  level  dependent  value;  given  here  at  E = EQ 


Table  V.  Comparison  of  <AE'>  (cm"  )  for  Different  Molecules 


O 

o 

o 

o 

O 

o 

to 

m 

o 

oo 

o 

VO 

00 

r>» 

f— 

CO 

o 

o 
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Appendix 

Vibrational  frequencies  and  parameters  for  RRKM  calculations  (ctn~b 
Nitromethane  molecule: 

3048(2),  2965,  1582,  1488,  1449,  1413,  1384,  1153,  1097,  921,  647,  599, 
476,  internal  rotation  (o  =  6). 

Activated  complex: 

3048(2),  2965,  1582,  1488,  1449,  1413,  1384,  647,  200(2),  120(2),  internal 
rotation  (a  =  6). 

E  =  54.4  kcal  mole”^ 
o 

F  =  1.3 


~r- 
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Fig.  1  First  order  law  plots  of  In  (a-x)  versus  ct  (sec) 

□  ,  T  =  870  K,  c  =  10'5  ;  A  =  T  =  1068  K,  c  =  10'2 


Fig.  2  Arrhenius  plots  of  log^  k  vs  1/T  for  the  m  =  5.0  and  m  =  14.4 
reactors. 

Fig.  3  Plots  of  Pc (m)  versus  T(K)  for  the  m  =  5.0  and  m  =  14.4  reactors, 
fl  are  the  experimental  points  for  the  m  =  5.0  reactor; 

C  are  the  experimental  points  for  the  m  =  14.4  reactor; 

-  end  -----  are  theoretical  plots  of  BE  model  and  FG  model, 

fitted  at  m  =  5.0,  for  the  m  =  14.4  reactor.  The  theoretical  curve 
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